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Recent results on muonic hydrogen [1] and the ones compiled by CODATA on ordinary hydrogen 
and ep-scattering [2] are 5a away from each other. Two reasons justify a further look at this 
subject: 1) One of the approximations used in [1] is not valid for muonic hydrogen. This amounts 
to a shift of the proton's radius by ~ 3 of the standard deviations of [1], in the "right" direction of 
data-reconciliation. In field-theory terms, the error is a mismatch of renormalization scales. Once 
corrected, the proton radius "runs", much as the QCD coupling "constant" does. 2) The result of 
[1] requires a choice of the "third Zemach moment" . Its published independent determination is 
based on an analysis with a p- value -the probability of obtaining data with equal or lesser agreement 
with the adopted (fit form-factor) hypothesis- of 3.92 x 10 -12 . In this sense, this quantity is not 
empirically known. Its value would regulate the level of "tension" between muonic- and ordinary- 
hydrogen results, currently at most ~ 4a. There is no tension between the results of [1] and the 
proton radius determined with help of the analyticity of its form factors. 

PACS numbers: 31.30.jr, 12.20.-m, 32.30.-r, 21. 10. Ft 



I. INTRODUCTION 

The results of a measurement by Pohl et al. [T| of the 
Lamb shift in muonic hydrogen and those compiled by 
CODATA on ordinary hydrogen and ep-scattering [5] are 
~ 5 a away from each other. The authors of £Q conclude 
"Our result implies that either the Rydberg constant has 
to be shifted by 2110 kHz/c (4.9 standard deviations), or 
the calculations of the QED effects in atomic hydrogen or 
muonic hydrogen atoms are insufficient." I discuss why 
the second option is part of the resolution of the apparent 
conundrum, but not all of it. 

It is intrepid [5] to use a model of the proton -in [T], 
a dipole form-factor- to challenge very well established 
physics -such as QED [TJ 0]. But this is not the only 
bone of contention: 

One of the approximations used in the theory of or- 
dinary or muonic hydrogen involves the lepton's wave 
function at the origin. The approximation is sufficiently 
good for the former atom, but not the latter. The re- 
quired correction can be rephrased by having an r p that 
runs, in the same sense as a s -the fine structure constant 
of QCD- does. The modification results in a ~ 3<r(/ii/) 
shift of the extracted central value of r p , in the direction 
of reducing the "tension" between experimental results. 
This correction depends on the model of the proton's 
charge distribution, but the model-dependence is a small 
correction to a moderate correction. These issues are 
discussed in detail in Sections III and IV. 

The current way to extract r p from ep scattering 
data involves an extrapolation to a momentum trans- 
fer q 2 = 0, the point from which (r p ) is inferred. This 
extrapolation covers a two-orders of magnitude larger 
hiatus than the one relevant to muonic hydrogen; the 



model-dependence is correspondingly larger. The extrap- 
olated object is a form-factor fit to data gathered above 
|q 2 | = 0(m 2 ), a domain where there is still "structure", 
relative to, e.g. a dipole form factor [5]. 

The extraction of r p from ep data has severe statis- 
tical problems, mentioned in the abstract and discussed 
in Section V. One way to reappraise this issue is to take 
new, very precise data [5], see also Section V. The diffi- 
culties associated with these analyses are shared by the 
measurement of the other relevant quantity: the "third 
Zemach moment" , as discussed in Section VI and VII. 

Sections VIII is a discussion of the experimental and 
theoretical results. Section IX contains my conclusions. 



II. THE ISSUE 

Former precise measurements of r p had two origins. 
One is mainly based on the theory [B] and observations 
[7] of hydrogen. The result, compiled in CODATA [2], is 

^(r2)(CODATA) = 0.8768 ± 0.0069 fm. (1) 

The second type of measurement is based on the theory 
and observations [3[9] of very low-energy electron-proton 
scattering. An analysis of the world data as of a few years 
ago yielded [5]: 

J{rf){ep) = 0.895 ± 0.018 fm (2) 

A recent ep-scattering experiment [5] results in: 

y^HAl) = 0.879(5) stat (4) syst (2) modcl (4) group fm, (3) 
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whose various subindices will be clarified anon. 

The proton's charge distribution, p p (r), is related to 
the non-relativistic limit of the electric form-factor, Ge, 
by the Fourier transformation 



G B (-q 2 )= / d\p p (v) 



1 - 



(4) 



which also serves to define (r 2 ) in proportion to the q 2 - 
derivative of the form factor at q 2 = 0. 

The most precise relevant measurement to date is that 



of the 2P*\= 2 



2S^ 2 1 Lamb shift in the pp atom, pQ: 



206.2949 ± 0.0032 meV. 



-\)\>. — >->■> <•->«.•>_ >m > ■ (5) 

In meV units for energy and fermi units for the radii, the 
predicted value [TU] is of the form 

L^[(rl),(rl) {2) ] = 
209.9779(49) 



5.2262 (r 2 ) 



0.00913 <r£) (2) (6) 



The first two coefficients are best estimates of many con- 
tributions while the third stems from the n = 2 value of 
an addend (6j E] 

AE 3 (n,l) = ^mt.S l0 {rl) (2) , (7) 



^m r Sio (r p ) (2 ) 
proportional to the third Zemach moment 
d 3 r 1 d?r 2 p{r 1 )p{r 2 )\r 1 - 



(8) 



For a specific model of p p (r) or its corresponding 
Ge{~ q 2 ), the two r-moments in Eq. |6| are related. For 
instance, for a dipole form factor 



A 2) 



(3675/256) [(r 2 )] : 



(9) 



while for a single pole [(?" 3 )(2)] 2 = (50/3) [(r 2 )] 3 . 

The authors of [T] use the dipole relation of Eq. ^ in 
Eq. ([6| to convert Eq. |5|, into an impressively accurate 



2 ) OH) = 0.84184 ± 0.00067 fm 



(10) 



The value of r p (pH) in Eq. (10) differs by ~ 3c(ep) 
from Eq. 5.0 cr(CODATA) from Eq. (I), and a bit 
more from Eq. ([3]). The standard deviations of these 
last three r p determinations are much bigger than the 
ones in Eq. ( |10[ ). Thus, they essentially determine the 
significance of the "distance" to the latter result. 



III. INSUFFICIENTLY-GOOD 
APPROXIMATIONS 

Let £ stand for e, p and let m r = me m p /(me + m p ) be 
the reduced mass. In an £p atom the dominant contribu- 
tion (99.45% of the total for I = p) to the coefficient of 
the (r 2 ) term in Eq. (|6|) is the familiar: 



- AE 



FS 

{n=2; 1=0} 
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<K) (ii) 



Recall that, in writing Eq. (11), the Fourier transform 
(V = — 47ra/q 2 ) of a Coulomb potential (V = —a/r) 
has been modified by the expression in the rhs of Eq. Q 
to obtain an additive term, oc 5(f), resulting in the "0" 
in the argument of the atom's wave function VF 

Even for I = p, the Bohr radius — l/(am r ) is 
orders of magnitude larger than r p , apparently justify- 
ing the consuetudinary approximation used in the last 
paragraph, which results in the ^(0) factor. But the 
precision of the measurement in Eq. ([5]) and its allegedly 
consequent Eq. ( 10 ) is so unprecedented, that the ap- 



proximation must be revisited, as I proceed to do. 

Consider a dipole form-factor Ge{— q 4 ) = m d/( m d + 
q 2 ) 2 , for which to 2 , = 12/(r 2 ). Repeat the analysis lead- 



ing to Eq. ( |ll[ ), this time without making the approxi- 
mation of Eq. Q . The result is 



" ZA - G {n=2; 1=0} 



1 + 3 m d a B + 4 m d 




(12) 



Naturally, the leading term coincides with Eq. ( 11 ). The 
first order correction to r p , estimated by entering the r 2 



value of Eq. (10) amounts to 0.42%. This may look tiny. 



But it increases the value of r p , extracted as in Eq. (10), 



by 2.7 a(pH). This modification of the central value of 
r p , though also insufficient by itself, is in the direction of 
reconciling the body of experimental results. 

It is also instructive to consider, for the nonce, a single- 
pole form-factor Ge(— q 2 ) = to 2 /(to 2 + q 2 ), for which 
2 _ 6/(r 2 ). The result is 



77J 
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(13) 



This correction amounts to 0.48%, or 3a(pH). Substi- 
tute to m for to 6 to conclude the obvious: for ordinary 
hydrogen and the precision of the corresponding obser- 
vations, the corrections of Eqs. ( 13|12 ) are negligible. 

We have learned that, at the level of accuracy of the 
pp experiment, the evaluation of the (r 2 ) term in Eq. rt5 



is not only delicate; it is also model-dependent. This is 
because of the inevitable extrapolation to q 2 = 0, where 
the radius is defined. We shall see that in the extraction 
of information from ep experiments, for which the ex- 
trapolation covers a two-orders of magnitude larger gap, 
the model-dependence is correspondingly larger. 



IV. A RUNNING {r 2 v ) 

The "atomic" subtleties discussed in the previous sec- 
tion are very familiar in QCD. To discuss the simplest 
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analogy, consider the total cross section for e + e~ anni- 
hilation into hadrons, above or in-between quark thresh- 
olds. It is of the form a(Q 2 ) oc (1/Q 2 )(1 + ajn + ...). 
For the approximation to be correct at all Q 2 , a s must 
"run", that is, be <5 2 -dependent in a specific way. 

In the next simplest example, the n-th moment of 
a (non-singlet) proton structure function -analogous 
to (r™)- if evaluated at two Q 2 values, differs by 
a multiplicative factor: to leading order, the ratio 
a s (Q 2 )/a s (Q' 2 ) to a specific anomalous dimension, d n . 



In a field theory, an expression like Eq. ( 11 ), containing 
a and an (r 2 ) referring to two different scales, would 
be a "mistmatch of renormalization points" . To correct 
it, one must evaluate ^ at the correct distance scale (as 
in the previous paragraph) or let (r 2 ) run. For a chosen 
form-factor this stat eme nt can be made precise, e.g. even 
the "5/-\/l2" in Eq. (12) has some meaning. 



The proton is not probed by the orbiting muon at r = 
0, or by momenta with equal weights in the range (0, oo) 
in the Fourier transform of 5(f). It is only probed by 
momenta ranging from |q| = 0(am r ) up to |q| ~ to<j/4, 
the proper "ultraviolet" scale. To use the mid expression 
in Eq. pTj ) at a consistent distance scale, r — 0, one may, 
for a dipole, define (fp)\( a m r ,md) = m d/^, keep the term 
| ^ (0)| 2 and substitute (r 2 ) by a running radius 



(^)l(o. 



oo) 



( r p) I (am r ,ra d ) 



1 + 5am r ,/(r 2 )| (Qmj . imd) /12 



(14) 



very reminiscent of the expression for a s in QCD. It is 



the lhs of Eq. ( 14 1 that is needed to extract the slope of 
the form factor at q 2 = 0, as in the rhs of Eq. Q. 

The dipole form factor is not foreign to QCD. The un- 
derstanding of the relatively high-Q 2 physics summarized 
by a dipole approximation, and the deviations thereof 
P21[T3] -as well as the related first measurement of Aqcd 
[12j- were discussed immediately after the discovery of 
QCD's asymptotic freedom [T4"] . 



THE EXTRACTION OF (rfj FROM ep 
SCATTERING DATA 



The Lyman-shift result quoted in Eq. JlO) is ^3a(ep) 
away from the ep-scattering result of Eq. (2l) . This is 
not a severe problem. A look at the data, reproduced 
in Fig. [I] on which the latter result is based, indicates 
that the problem if even less severe. What is shown in 
the figure are data available in 2003, normalized to a 5- 
parameter continued-fraction expansion of Ge(— q 2 ) [Hj- 
The fit's result is x 2 / n doi — 1-652, or, more explicitly, 
X 2 — 512 for n do f — 310 degrees of freedom. 

It may be useful to recall that the p-value of a data- 
set relative to a given assumption or fit -in this case 
the specified continued fraction- is the probability of ob- 
taining data at least as incompatible with the hypothesis 
as the data actually observed. Let /(x 2 ,ndof) be the x 2 
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FIG. 1: Low-|q| data, compiled and analysed in [51 111]. 



probability distribution function. Let T(a, b) \T(b)] be the 
incomplete [ordinary] gamma function. Then 



P(X i«dof) = / f(z,n dof )dz = 



r(n dof /2, X 2 /2) 
r(n dof /2) 



(15) 



and p(512, 310) ~ 3.92 x 10~ 12 , i.e. the quality of the fit in 
[5] is not "quite good" . It is possible [5] to reduce this be- 
hemoth disagreement by adding quadratically 3% to the 
Stanford error bars (to obtain p(370, 310) ~0. 011), or by 
a norm change of 1% of [TS]... [which] would decrease x 2 
by 60 (resulting in p(452, 310) ~2.38 x 10" 7 ). Modifying 
the data is not necessarily a universally accepted proce- 
dure, or so would the corrected experimentalists opine. 

It is also possible to draw sensible-looking curves 
through the data that, in their slope at q 2 = 0, differ 
from a straight horizontal line in Fig. [T] by one or more 
of the ct's in Eq. pj). The fact that the data points are 
very scattered is an unavoidable problem. One way to 
reconsider the issue is to take new and very precise data. 

Such data exist [S] and are partially reproduced in 
Fig. [2j The paper contains many relevant commentaries. 
One of them is: "The structure at small Q 2 seen in Ge 
and Gm corresponds to the scale of the pion of about 



Q 



0.02 (GeV/c ) and may be indicative of the 



influence of the pion clowd." The most apposite remarks 
in [S] concern the extraction of the results: 

Two types of "flexible" models are considered in [S]: 
fits to polynomials and spline fits. The Ge results are 

r 2 ) (spline) = 0.875(5) st at(4) S y S t(2) mo doi fm, 

r 2 )(polynom) = 0.883(5) sta t(5) sys t(3) m odci fm (16) 

"Despite detailed studies the cause of the difference could 
not be found. Therefore, we give as the final result the av- 
erage of the two values with an additional uncertainty of 
half of the difference" [5j : the outcome quoted in Eq. (J3j) . 
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FIG. 2: The lowest-energy Mainz data [5|. The dark trumpet 
encompasses estimated theoretical and systematic uncertain- 
ties and ±1ct statistical errors. Barred points are results of 
other experiments. The dotted lines are by-eye fits to data in 
the domain AQ 2 = (~ 0.06, 0.2) (GeV/c) 2 , for each of the sets 
of data. Their arrowy continuations are meant to illustrate 
how difficult it may be to extract from such an extrapolation 
the slope at the lowest-Q 2 measured point. 



THE EXTRACTION OF <r£) (2) FROM ep 
SCATTERING DATA 



The result for the third Zemach moment is [TT] 

[(^}(2)] 1/3 = (1.394 ± 0.022) fm, (17) 

based on the data in Fig. [T] Based on the same data, the 
result for (r 2 ) is that of Eq. ((2j). To discuss the "third" 
moment, it is useful to write it in an alternative form: 



^(q 2 ) = 



o 

48 



d|q| /(q 2 ) 
G|(q 2 )-1 
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(18) 



Notice that /(q 2 ) tends to a constant as q 2 — > 0. 

The usual dipole form factor (m 2 , = 0.71 GeV 2 ) is a 
good-enough approximation for the forthcoming discus- 
sion. The shape of /(q 2 ) is shown in Fig. [3j Normalized 
to the total integral, the fractions of the integral in vari- 
ous relevant intervals Aq = (o, b) fm -1 are the following: 



(0,0.4)^0.09; (0.6, 1.6) ->• 0.20 ; 
(1.6,4) 0.29; (4,oo)^0.38. 



(19) 



Whether the fits' uncertainties are thus correctly esti- 
mated is debatable, but this is not the main point. 

The crux of the matter is that the procedure in [5] illus- 
trates how, even for sets of "flexible" fits, the result is sig- 
nificantly set-dependent. The reason is simple: two ana- 
lytic (or piece- wise analytic) functions arbitrarily close to 
each other in a given interval, say AQ 2 = (<5 2 nin , oo), can 
be arbitrarily different in their continuation to another 
interval, such as AQ 2 — (0, Q 2 ^). 



The data itself could be used to study the model- 
dependence of the extracted value of (r 2 ). Suppose that 
one fits the data in the interval AQ 2 = (~ 0.06,0.2) 
GeV 2 of Fig. [2] and extrapolates to the lowest-q 2 point 
at which there is still data and the Ge slope is measured. 
This is analogous to extrapolating to q 2 = 0, except in 
that the answer is known. A look at Fig. [2] suffices to con- 
clude that the result is likely to be significantly wrong. 

The "less flexible" models used to analyse the Mainz 
data have x 2 / n dof ~ 1.16 to 1.29 for n dof w 1400 [5]. 
The corresponding p-values range from 2.69 x 10~ 5 to 
9.07 x 10~ 13 . The most flexible ones have x 2 /n do f ~ 1.14, 
or p = 1.88 x 10 -4 , that is, they are far from being flexible 
enough to describe the data. There seems to be a general 
tendency to forget that the quality of a fit is a function 
of two variables, not of their ratio, and that for large 
fixed ndof the dependence of the fit's quality on \ 2 is an 



inordinately sharp function around \ 
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to plot Eq. ( 15 ) to convince oneself. 



FIG. 3: The third Zemach integrand of Eq. (18 1, illustrated 
for the usual dipole and shown in the same domain as that of 
Fig. 111. The shaded band is at Aq = 1.1 ±0.5 fm -1 . Data in 
the dashed extrapolation do not exist. The contributions to 
the moment from various intervals are discussed in the text. 

It is stated in |llj that "Sensitivity studies have shown 
that the main contribution to the integral comes from 
the region Aq = (0.6, 1.6) fm -1 where the data base for 
electron-proton scattering is very good." But the con- 
tribution is significantly larger in the range Aq x fm = 
(1.6,4) where the data are particularly bad, see Fig. I] 
The contribution if the range (4, 00) is even larger. The 
contribution in the range (0, ~ 0.4), where there are sim- 
ply no data, is not at all negligible. In other words, given 
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the results in Eq. (19 1, the quoted statement must refer 



to the error estimate, not to the central value. And the 

r 3 ) 
pi 



basis for the deduced central value on (r„)( 2 ) is still the 



fit in Fig. [T] whose p- value I have quoted. 

It is concluded in another study [TC] of /(q 2 ) that "a 
large third Zemach moment can only occur if (r^) is also 
large". This is unobjectionable, though "large" means 
relative to the expectation from a dipole form-factor. 



VII. HOW TO EXTRACT 



MOMENTS FROM 



ep SCATTERING DATA? 



We have seen that the values and errors obtained in 



the literature for 



and 



( r p){2) are n °t credible. 



Since 



the various moments are highly correlated, another perti- 
nent question is: how to draw, in the ((r 2 ), (t" 3 )( 2 )) plane, 
trustable, model-independent contour plots of given sig- 
nificance? A procedure might be the following: 

(1) Normalize the data to a fixed model, as in Figs. [l]2] 

(2) Study modifications, relative to the model, with 
a complete set of orthogonal functions, e.g. a discrete 
Fourier basis for the complete data interval, a function 
of a variable, such as logq 2 , chosen to emphasize the 
most relevant, low-q 2 , domain. 

(3) Let the results fix the needed flexibility, i.e. cut the 
Fourier series at the term for which \ 2 *S ^dof (p *S 0.5). 

(4) Sidestep an extrapolation to q — > 0, which is un- 
avoidably problematic. That is, use the data only where 
they exist. For this, one would have to Fourier transform 
Gg(q) into p(r) and study its moments. This is probably 
the only way of facing their unavoidable correlations. 

(5) Show the correlated results as ((r 2 ), (t" 3 )( 2 )) contour 
plots for fixed acceptable values of p. 

Such a procedure is very different from the usual. It 
may well lead to significantly different conclusions. Doing 
this analysis -in contrast to the simpler choice of verbally 
discussing it- is well beyond the scope of this paper. 



VIII. DISCUSSION 

The result of Eq. ([I]) is not only based on measure- 
ments including ordinary-hydrogen levels, but also on the 
ep-scattering result of Eq. ^2y, discussed in Section V. 
The elimination of this input from the COD ATA fit to 
78 (more or less) fundamental constants results in [2J: 



r 2)( ep ) = 0.8737 ± 0.0075 fm 



(20) 



This result is shown as the shaded band in Fig. |4j dis- 
playing the rms radius versus the cubic root of the third 
Zemach moment. To facilitate the coming discussion, 
I have added the — 2er and — 3cr lines corresponding to 
Eq. (20 1. Also shown in the figure are the two results, 
Eqs. (161, of the Mainz experiment 0. 



The lowest (dashed) line in Fig. 4] is Eq. (10), from 
the \xli Lyman shift [T]. The continuous straight line 
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FIG. 4: Atomic data for the rms radius versus the third 
Zemach moment. The curve is the correlation between 

the rms radius and the third Zemach moment implied by the 
observation of Eq. Q and the theoretical result of Eq. (21 1. 
Also shown are the polynomial and spline results of 5i (meant 
to be read as horizontal bands) and the (rp}( 2 ) of 



above the previous one takes into account the renormal- 
ization correction of Eq. (13). Make the same correction 
in Eq. ([6| to obtain: 



209.9779(49) - 5.20123 (r 2 ) + 0.00913 (r 3 ) (2) (21) 

and equate it to the observed value of Eq. ([5| to obtain 
the "/iif" correlation, the continuous curve in Fig. [4] 
This correlation -and not a figure for r p — is the outcome 
of the theoretical analysis of the measurement pQ. 

If the value of the third Zemach moment was that of 
Eq. (17), indicated by an arrow in Fig. [4j the muonic- 
and ordinary-hydrogen results would be more than 3<r 
away. If, instead, ((j" 3 )^)) 1 / 3 ~ 2.1 fm the tension would 
diminish to the 3cr-level, marked by a circle in Fig. [4] (the 
increase of the moment is more severe than it seems to 
be, since the observable is not its cubic root). 

The standard deviations of the previous paragraph are 
the ones pertinent to a normal data distribution, for 
which ±1, ±2 and ±3er correspond to coverage proba- 
bilities of 68.27%, 95.45% and 99.73%. But the data of 
CODATA are not normally distributed, meaning that the 
bands of the same fixed probability are not the ones in 
Fig. [4] and that the conclusions of the previous paragraph 
should be correspondingly weakened. 

More precisely, 9 out of 135 input data in [2] -related 
to the Watt balance, the lattice spacing in various Sili- 
con crystals, the molar volume of the same element and 



G 



the quotient of Plank's constant to the neutron mass) 
have had their uncertainties increased by a multiplica- 
tive factor 1.5 |1] . This choice helps in obtaining a fairly 
satisfactory overall p-value, p = 0.221, but it describes 
a hypothetical set of experiments, not the actual one. 
Moreover, the question of the individual p-values of the 
experiments is not reexamined in [2]. We have seen ex- 
amples of how misleading this omission may be. 
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FIG. 5: The CODATA result of Eq. |20f , the fj,H correlation 
also shown in Fig. [4j and recent "theory-driven" results for 
(r-p) (meant to be horizontal bands). Bl, B2 cite [17] . HP1 to 
HP4 cite [T5] and WLTY quotes P2]. 



A solidly-motivated approach to the extraction of (r^) 
from the ep scattering data is based on the use of the 
analytical properties of the nucleon's form factors. Two 
recent examples are \F?\ [18] (it would be helpful to know 
the p- values of these fits). Their r p outcomes are shown 
in Fig. [5j The result of a first-principle lattice calculation 
is also shown. The spread of the theoretical results 
may be indicative of the difficulty of reaching a consensus. 
Should the spread reflect a level of uncertainty, there is 
no "tension" between theory and observations. 



this is because the wave function is not only probed at 
the origin, but up to distances at which details of the 
proton charge distribution (beyond its mean square ra- 
dius) become relevant. Thus the model-dependence of 
the correction, see Eqs. ( p|l2| ). 

Translated into momentum space, the previous para- 
graph becomes very familiar. A wave function at the 
"origin", a 6(0), corresponds to a uniform sampling of all 
momentum transfers, something that an atom can hardly 
provide. The correction, a proper match of renormalisa- 
tion scales, can be phrased as a proton's running radius. 
Not a surprising result: in a field theory all measured 
quantities are scale-dependent. 

The above correction to the analysis in IF reduces its 
disagreement with the (r^) CODATA result from ~ 5<r 
to at most ~ 4cr. The "at most" is crucial, for the 
conclusion depends on the value adopted for the third 
Zcemach moment. We have seen that its extraction from 
ep scattering data is most questionable. In Fig. [4] one 
can see that, even after the corrections I have discussed, 
(( r p)(2)) 1 ^ 3 ~ 3 fm would be required to have ordinary 
and muonic hydrogen precisely agree. Even if a dipole 
form factor is only a very vague description of the data, 
such a value feels unexpectedly large, part of the argu- 
ment in |16j . 

The Lamb shift measurement provides a correlation 
between the two relevant moments: the narrow curved 
domain labeled "//if" in Fig. [3] It would be very help- 
ful to extract the correlation dictated by ep data, to be 
added as confidence-level contours to the figure, to de- 
cide -with confidence- what the empirical conclusion is. 
It may well be that the /iH and ep correlations have a 
sufficient overlap for the question of data incompatibil- 
ity to be moot. After all, also for ep scattering, the two 
quantities plotted in Fig. [4] are obviously strongly corre- 
lated, a fact that has been totally ignored. 

Similar inferences can be extracted from the compar- 
ison of theory and data summarized in Fig. [5j Cur- 
rently none of these "theory-driven" results are available 
in the form of two-dimensional ((r^), (fp)(2)) confidence- 
level plots. Even barring other putative limitations of 
current theory or experimental analyses, the most ex- 
treme views on the subject at hand [31 2] seem to have 
been largely exaggerated. 



IX. CONCLUSIONS 

We have seen that, to the precision required to analyse 
the recent muonic-hydrogen results, the vintage "wave- 
function at the origin" expression for the Lamb shift, 
Eq. (11), is insufficiently precise. In configuration space 
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